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Abstract 



> . 

VO . The Dyson- Schwinger equation for the quark self energy is solved in rain- 

bow approximation using an infrared (IR) vanishing gluon propagator that 
introduces an IR mass scale h. There exists a h dependent critical coupling 
^^ indicating the spontaneous breakdown of chiral symmetry. If one chooses 

Q>^ ■ realistic QCD coupling constants the strength and the scale of spontaneous 

(— I , chiral symmetry breaking decouple from the IR scale for small b while for 

Qh! large h no dynamical chiral symmetry breaking occurs. At timelike momenta 

Q^' the quark propagator possesses a pole, at least for a large range of the param- 

eter b. Therefore it is suggestive that quarks are not confined in this model 
for all values of b. Furthermore, we argue that the quark propagator is ana- 
lytic within the whole complex momentum plane except on the timelike axis. 
Hence the naive Wick rotation is allowed. 
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I. INTRODUCTION 

New accelerators like CEBAF, MAMI-B and COSY will investigate liadron observables 
at a scale intermediate to the low-energy region where hadron phenomena mainly reflect 
the underlying principles of chiral symmetry and its dynamical breakdown and to the high 
momentum region where the "strong" interaction has the appearance of being a perturba- 
tion on free-moving quarks and gluons. A reasonable theoretical description of intermediate 
energy physics therefore has to satisfy at least the following requirements: Formulated in 
terms of quarks and gluons it has to account for a mechanism of dynamical chiral symmetry 
breaking (DCSB) and identify the related Goldstone bosons in the measured spectrum of 
quark bound states. In addition, all the correlation functions have to transform properly 
under the renormalization group of Quantum Chromo Dynamics (QCD), i.e. at large mo- 
menta one should recover the correct anomalous dimensions. Only if these requirements are 
fulfilled the theoretical framework is fast in those grounds of the theory of strong interaction 
that are well established. Furthermore, it is desirable to formulate a microscopic picture of 
the cause of confinement. 

In this context growing interest has been focused on the Euclidean Dyson-Schwinger 
equations (DSE's) as a tool for developping and applying non-perturbative methods in quan- 
tum field theories, a recent review on this subject is given in Ref. |l[]. Realistic models of 
QCD derived through truncation of the DSE tower are applicable in the whole Euclidean 
momentum space and, among others, have the following features: Quark condensation gives 
rise to DCSB, the chiral Goldstone bosons are identified with the lowest mass pseudoscalar 
quark- ant iquark bound states, perturbative QCD is matched at large momenta, and a study 
of the analytic structure of the quark propagator provides some insight into the mechanism 
of confinement. 

The DSE's include the Bethe-Salpeter equations (BSE's) for bound states, especially the 
ones describing physical mesons ^. The BSE's are coupled to the DSE's for the quark 
and gluon propagators, i. e. the BSE's need as input at least the full renormalized quark and 
gluon propagators which — for calculating the meson mass spectrum, for instance — have to 
be known at complex Euclidean momenta far away from the real axis. In principle, the quark 
and gluon propagators can be determined from their DSE's. As these are coupled to the 
infinite tower of DSE's one is not able to solve the equations exactly. Possible approximation 
schemes are based on the truncation of the DSE tower: Some n-point functions like the 
gluon propagator and the quark-gluon vertex are parametrized choosing Ansdtze such that 
certain requirements are fulfilled. Important constraints are imposed by the discrete and 
continuous symmetries of QCD most of which are formulated in Ward or Slavnov- Taylor 
identities. Additional input may be taken from hadron phenomenology. 

In this paper we investigate the analytic structure of the quark propagator using the 
Euclidean DSE for the quark self energy in rainbow approximation where the full quark- 
gluon vertex is replaced by the perturbative one. The main ingredient for our study is a 
color diagonal gluon propagator whose transverse part vanishes at zero four momentum: 

^"(*^) - F^ ■ (1) 

Herein k denotes the (spacelike) gluon momentum and b is an energy scale to be determined. 
This form of the gluon propagator is motivated by very distinct considerations. One is based 



on the complete elimination of Gribov copies within a lattice formulation of non-abelian 
gauge theories |^: Landau gauge is fixed uniquely with the help of a 'thermodynamic 
parameter'. The resulting functional integral is dominated by field configurations on the 
Gribov horizon. This implies that the gluon propagator vanishes for small momenta k 
as k"^. In a couple of different studies a gluon propagator of the form (^ is obtained by 
considering a generalized perturbative expansion of the seven superficially divergent proper 
vertices of QCD which are allowed to depend non-analytically on the coupling constant [Q. 
Furthermore, the Field Strength Approach (FSA) to QCD supports the form (|l]) for being 
the dominant part of the infrared quark-quark interaction |^ ; however, the color structure of 
the FSA gluon propagator is not diagonal. At last, some recent lattice calculations obtained 
a gluon propagator which allow a fit of the form (|l]) (but which do not rule out a fit using 
a standard massive particle propagator) [^]. 

As the gluon propagator (|l|) has no poles for timelike momenta it may be argued that it 
describes confined gluons. Instead, the gluon propagator has poles on the imaginary axis, 
k"^ = ^ib"^, and therefore the gluon 'decays' after some time r ~ 1/b. This can be interpreted 
as the gluons being screened. On the other hand, the gluon propagator (|T]) is smaller than 
the perturbative one for all momenta: k'^ /{k^ + h^) < l/k"^. This has drastic consequences 
for DCSB as we will show: For the values of b that are suggested by the lattice studies of 
Ref. 0] chiral symmetry will be realized in the Nambu-Goldstone mode for unrealistic large 
coupling constants only. We also find that quarks are not confined in this model, at least 
for a large range of model parameters b. 

A similar study was undertaken by Hawes et al. [0]. In this work the gluon propagator 
(|I|) and a non-trivial quark-gluon vertex has been used to calculate the quark propagator 
in the spacelike region. Their results are similar to ours, especially they find unconfined 
quarks. In contrast to our work where the analytic continuation to timelike momenta is 
performed explicitly and a pole signaling free massive quasi quarks as asymptotic states is 
found, Hawes et al. use a test for confinement employing a Fourier transformation of the 
quark Schwinger functions. The large time behavior of these Fourier transforms indicate 
asymptotic states of massive deconfined quarks. 

In Sec. 2 we present our model DSE for the quark self energy. We also discuss the Ansatz 
Eq. (0) for the gluon propagator in more detail. In Sec. 3 we present the numerical solution 
for the quark self energy at real Euclidean momenta and discuss the dynamical breaking of 
chiral symmetry. In Sec. 4 we describe the analytic continuation of the quark self energy with 
special emphasis on its pole structure and the Wick rotation from Euclidean to Minkowski 
space. In Sec. 5 we conclude. 

II. DYSON-SCHWINGER EQUATION FOR THE QUARK SELF ENERGY 

Our starting point is the DSE for the QCD quark two-point function in Euclidean momen- 
tum space which can be written in terms of the quark self energy S(p). After renormalization 
it reads 0: 

S(p) = (1 - Z^h ■ V + {Z,Z^ - l)m + Z, Cpg^ J ^ D.^ip - q) T,{p, q) S{q) 7. (2) 



with Cf = (^c'~l)/(2^c) = 4/3 for three colors and g denoting the QCD couphng constant. 
m is the current quark mass. Zi(/i^,A^), Z2(/i^,A^) and Zm(/i^,A^) are renormahzation 
constants depending on a regularization parameter A and getting fixed at the renormahzation 
point p^ = /i^. 

The three Green's functions that appear on the r.h.s. of Eq. (|^) are 

1. the quark propagator, 

S{q) = 1 ^ 1 ^ Z{q') 

7 ■ g — m — S(g) 7 ■ qA{q'^) — B{q'^) 'J ■ q — M{q'^) 

with B and A denoting the scalar and vector parts of the quark self energy S, respectively, 
and Z and M the field renormahzation function and the mass function, respectively; 

2. the gluon propagator 

D,Ak) = [s,^ - ^) Drie) + e^ (4) 

where only the transverse piece, Dt, is dressed due to vacuum polarization and 

3. the quark-gluon vertex r^(p, g). 

Of course, the gauge parameter ^ appearing in the gluon propagator (^) is a function of 
the renormahzation point /i^, i.e. any redefinition of the renormahzation point leads to a 
new gauge. The only gauge which is not affected by renormahzation is the fix point ^ = 0, 
i.e. Landau gauge P,p!0|. In our studies we therefore concentrate on this case. 

For studying Eq. (Q) one has to put in the full gluon propagator and the full quark-gluon 
vertex; an exact investigation of the quark DSE (|^) therefore requires solving their DSE's, 
too. To do so the gluon 3- and 4-point functions, the ghost self energy, the quark-quark 
scattering kernel, and so on, have to be known. Each of them again fulfills its own DSE. We 
do not address the problem of solving all these DSE's simultaneously but instead truncate 
the infinite tower of DSE's by making Ansdtze for both the gluon propagator as well as the 
quark-gluon vertex. 

Neglecting ghost contributions to the gluon vacuum polarization the transverse part of 
the gluon propagator (|), Dxik"^), is closely related to the running coupling constant of 
QCD, a(P): 

/i,,(,=) = «. (5) 

In addition, 

Zi(Auv,/i') = Z2(Auv,/x') (6) 

in Landau gauge; this is a result of the Ward identity that connects the quark-gluon vertex 
to the inverse of the quark propagator. 

Taking into account one-loop corrections to the QCD perturbation expansion one possible 
smooth interpolation between IR and UV behavior is the following Ansatz for the gluon 
propagator: 

g^ D^e) = Df{e) , , i,.. , (7) 
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with Aqcd the QCD scale parameter and A the couphng strength which is related to the 
anomalous dimension of the quark mass, dM = 12/(33 — 2iV^) [|l^ , |T5|] : A = iii'^dM- However, 
in our numerical studies we will treat A as an adjustable parameter. Furthermore, we set 
r = e to ensure that the whole strength of the infrared (IR) quark-quark interaction is 
carried by the function D^{k^). In the following we will call the second factor on the r.h.s. 
of Eq. (^ 'ultraviolet (UV) improvement'. Note that in Eq. (|^) the four-momentum k is 
assumed to be spacelike. We like to mention that its continuation to large complex momenta 
is by no means fixed by the behavior in the asymptotic spacelike region. For example one 
could replace 

irdM 2'KdM , . 



In (r + P/A^cd) In (r^ + ^VA^cd 



without altering the spacelike UV behavior. Indeed, both expressions are quite similar even 
for small spacelike four-momenta. On the other hand, their analytic properties are very 
different. While the function on the left hand side has its cut on the timelike axis, k G ilR, 
the one on the right hand side has cuts for k G y/SlR. 

Of course, one can use a two- or a three-loop order running coupling of QCD instead of 
the UV-improvement in Eq. (^. Nevertheless, even these cannot provide information about 
the structure of the gluon propagator at small spacelike momenta (< 2 GeV). Hence, for 
proceeding a parametrization of D^ is unavoidable. 

On the parametrization of the IR part of the gluon propagator a great deal of work has 
been done during the last years. Several parametrizations have been investigated and have 
supported the opinion that quark confinement is closely related to infrared slavery or at least 
to a strong finite quark-quark interaction in the infrared momentum regime. Nevertheless, as 
remarked in the Introduction, there exist very distinct studies 0, 0] (using Landau gauge) 
and that propose 

Ofik^) - pf^ (9) 

with a dynamical mass scale h assumed to be proportional to Aqcd- In addition, Refs. [^,|^ 



even argue that h k, Aqcd- This propagator has no Kallen-Lehmann representation |T^ 
indicating that the propagation of a free quasi gluon is not possible. As the strength of 
the quark-quark interaction vanishes at small momenta (corresponding to large distances in 
coordinate space) the solution of the quark DSE @j with the gluon propagator (P) as input 
would exhibit a new mechanism for quark confinement if any. 

For simplicity we use the rainbow approximation, i.e. we formally replace the dressed by 
the perturbative quark-gluon vertex in Eq. (|^): 

r^(p,g)=7M- (10) 

In this approximation the solutions of the Bethe-Salpeter equations (BSE's) for color-singlet 
quark- ant iquark bound states show the correct low energy behavior as formulated in cur- 
rent algebra and partial conserved axial current (PCAC) relations. For instance, the quark 
bound states carrying the quantum numbers of the pseudoscalar mesons can be identified 



with the Goldstone bosons related to the dynamical breakdown of the SJjR^Nf) x SUiiNf) 
chiral symmetry of massless QCD |11]; in case of massive QCD the pion mass is deter- 



mined by a generalized Gell-Mann-Oakes-Renner relation. The rainbow and the separable 
approximation [l^ are the only known approximations that imply the foregoing consistent 
formulation of quark-DSE and BSE for mesons. Furthermore, the rainbow approximation 
(p!0|) is equivalent to the stationary phase approximation of the global color-symmetry model 
(GCM) []T^, an abelianized version of QCD with quark fields being transformed to collective 
fields composed of quarks and anti-quarks and carrying the quantum numbers of mesonic 
fields. 

After inserting Eqs. (|^) and (p!OD (together with or without Eq. (0)) into Eq. (0), using 
Eq. (^ and fixing the renormalization point, /i^, such that vector and scalar part of the 
quark self energy obey the renormalization conditions 



B{p' 



/i 



a 



(11) 
(12) 



with constants a and (3 being fixed later, we get two integral equations which couple the 
functions A and B that now depend on the momentum s = p"^, the renormalization point 
//^, and the regularization parameter A^y-. 



A{s,fi^;A 



uv) 



a 



with 



B{s, fi^;Aly) =(3 + a 



JA[A,M]{s,f^^;Aly) 



-JA[A,M]{s,fi^;Aly) 

J^[A,M](/i2,/i2;A2jv) ' 

Jb[A, M](s, /i^; Aly) - Jb[A, M]{fj.^, /i2; A 
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(13) 
(14) 
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dxlA^x; r, s) H(x; r, s) 



JB[A,M]{s,f^'';A^ 



A r^l 



drM{r,n'^]A 
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67r3 Jo A{r, fi^; A^^^) [r + M^r, /i^; A^^^)] 



dxlsix; r, s) S(x; r, s) 



(15) 



(16) 



where 



without 

with UV improvement, Eq. (7) . 



log (r + r + s - 2xy/¥^/Al^cj^ 

The functions I a and Is are defined in the appendix, Eqs. (|A1| ) and (^), where also the 
analytic expressions for the angle integrals (x-integrals) in case of S = 1 are given (Eqs. (|A3|) 
and ( |A4|) ). The mass function in Eqs. ([T3|) and (0) reads: 



M(s,/i";A^v 



E(g,/i^;A^v) 
A{s,fi^;Aly) 

p[l- JAJfi^ fi^ A^v)] + a [JbJs, /i^ A^v) - Jsifi'', fi^; A^y)] 
a[l- JA(s,/i2,A^v)] 



(18) 



To make contact with perturbative QCD one has to formulate an appropiate renormal- 
ization condition ensuring the asymptotic freedom of quarks. It is convenient to choose the 
renormahzation point as large as possible and set 

a = 1 , (19) 

(3 = m (20) 

at this point. The largest available momentum is s = A^v! therefore we set: /i^ = A^y if 
not stated otherwise. 

Numerically, it turns out that Jyi[y4, M](/i^,/i^; Ayy) is very small for /i^ being large 

,2 



enough. Actually, in the limit /i^ ^ cxd it vanishes exactly. Hence, solving the DSE (|T3D 
and (0) for the renormalized quark self energy is (within numerical accuracy) identical to 
solving the equations 

A{s, /.2; kly) = 1 - Ja[A, M] (s, /i^; A^y) , (21) 

5(s,/x2;A^y) = m(/i2) + j^[A,M](s,/i2;A^y) , (22) 

when imposing the renormahzation condition 

,./ 2 a2 ^ _ mA(/i^ /i^; Agjy) + B{s, /i^; A^y) - g(^^ /i^; A^y) 

m is a running current mass 

m{^^^)=m-JB[A,M]{^l\^^^■Kly). (24) 

We like to remark that in case of the DSE system (^,^) this condition is suf- 
ficient because the approximation J^(/i^, /x^; Auv) =0 sets the renormahzation constant 
Z2(/i^;Auv) = 1. For the system ([13| , p!^) Z2 has to be fixed dynamically and hence two 
constraints, Eqs. (|TT|) and (|12D , are needed. 

Asymptotic behavior of the quark self energy 

At large spacelike momenta s the IR mass scale h in the gluon propagator (^ can be neglected 
and the integral equations (|1^) and (|1^ with S = 1 [i.e. no UV improvement) can be 



converted into differential equations |2^. As for large s the integral Ja[v4, M](s,/i^, A^y) 
becomes very small the vector part of the quark self energy is approximatly equal to 1 and 
hence Z2 ~ 1. The scalar part B is then determined by the non-linear differential equation 

(fB ^dB \ B ^ . . 



with the boundary condition formulated in Eqs. ( [I^J20| ). If there exists a bounded solution 
one can drop the term B^ in the denominator and one gets the asymptotic solution 



Bis) s large a 

— r4 ^ , <; cos 



In [s//i^j + (p 



(26) 



with K = J\/\c — 1- Ac = vr^ ^ 9.87 is (up to a trivial color factor Cp^ = 3/4) the critical 
coupling of Quantum Electro Dynamics (QED). Implying the UV boundary condition (pOD 
one gets 



(y9 = arccos — . (27j 

a 

The parameter a is uniquely fixed by the IR boundary condition i. e. by the solution of the 
integral equation (|Iip which determines the slope of B at the renormalization point /i^. 

In case of the UV improved model as defined by the gluon propagator (|^) the asymptotic 
behavior of the mass function M{s) differs qualitatively from Eq. (|26|) . An analysis based 
on the renormalization group yields [p^,p]] : 



s large A , , 1 



log (/iVA 
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M{s) > -o(gg).^T . ) . La/4.^ - (28) 



log (s/A 
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Dynamical Chiral Symmetry Breaking 

Compared to the hadronic mass scale Aqcd the up and down quark current masses m are 
negligible. In the limit of zero current masses, m = 0, the QCD Lagrangian is invariant 
under chiral transformations. The Wigner-Weyl realization of the vacuum state corresponds 
to the trivial solution of Eqs. (p!3| ) and (|T^: 

5triviai(s,/i',A^v) = 0. (29) 

In QED in Landau gauge this trivial solution maximizes the CJT action |T^ whose stationary 
phase condition is identical to the quark DSE (0) [0,^. This statement holds also true 



for abelianized QCD as long a.s A ^ 1. Hence, a vacuum configuration with B ^ will 
be dynamically favored ii A ^ 1. This may be tested by evaluating the CJT action at the 
stationary points. 

A nontrivial solution of Eqs. ( ]I3| ) and (|T^ (in the case m = 0) indicates the dynamical 
breaking of chiral symmetry, i.e. the corresponding vacuum state is realized in the Nambu- 
Goldstone mode. From the definition of the quark condensate 

it is obvious that a non- vanishing scalar part of the self energy, B > 0, leads necessarily 
to a non- vanishing condensate {qq) ^-2 < 0. The scalar part of the quark self energy at zero 
momentum, -B(O), or, correspondingly, the mass value M(0) is therefore an order parameter 
of DCSB. 

Confinement and the analytic continuation of the quark self energy 



The analytic continuation of Eqs. (13) and (|T4D to complex momenta is very important in 

the studies of DSE's for the following reasons: 

1. A quark propagator whose mass function obeys the relation 

sm + M\sm)=^ (31) 

and whose renormalization function Z is non-vanishing at s = sm < describes a massive 
particle which should be detectable in some appropriate experiment. Hence, one way to 
obtain information about quark confinement is the examination of the pole structure of the 
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quark propagator after continuing real Euclidean to imaginary momenta, i.e. s to —s: If 
there does not exist any momentum sm for which Eq. (31) is fulfilled, there is no Kallen- 



Lehmann representation for the quark propagator and the described quasi-quark is confined. 

2. In order to solve the quark bound state equations {i.e. the BSE's for mesons or the 
Fadde'ev equations for baryons) information about the quark self energy in a large domain 
of the complex momentum plane is needed. While external momenta {e.g. the momentum 
of a meson close to its mass-shell) can be timelike the loop momentum occurring in the 
BSE is spacelike. Hence, there appear combinations of spacelike and timelike momenta as 
arguments of the quark self energy functions A and B; the analytic continuation of these 
functions to complex momenta is therefore unavoidable if one wants to describe hadrons 
others than pions. 

3. Studying the analytic properties of the quark self energy provides some insight into the 
connection of the Euclidean formulation of QCD and QCD in Minkowski space. Recently 
it has been shown [^ that in models containing both gluon and quark confinement the 
Wick rotation cannot be performed naively indicating thereby the close relation between 
confinement and singularities or branch cuts appearing in the quark propagator at some 
complex Euclidean momenta. 

III. QUARK SELF ENERGY FOR SPACELIKE MOMENTA 

We have solved Eqs. (0) and (|1^) by iteration with fixed cutoff Auy. Different 
parametrizations and grid densities of the discretized momentum interval s G [0, A^v] have 
been used. All the results we report herein do not depend on changes of the grid parametriza- 
tion or further increasements of the grid density. 

To check the accuracy of the numerical procedure we first study the asymptotic behavior 



of the DSE-model without UV- improvement, i.e. we set S = 1 in Eqs. (|T3D and (jT^). At large 
renormalization points we fix the mass function arbitrarily and fit a tail of the form (^Hf) 
in the momentum region (6^ <^) /i^ < s < A^y allowing all parameters a, (p and Ac to get 
varied. 

In Fig. 1 the function B{s)/B{0) is shown for A = 16, /i = 2^^6, Ayy = 2^^b and 
m/B{0) = 0.0153 together with an asymptotic fit (^) with a = 0.06858 and 

ip = 1.3461 , Ac = 9.869 . (32) 

A comparison with the calculated values Ac = tt^ ~ 9.870 and ip = arccos(m/a) = 1.3461 
shows high accuracy for the numerical solution of Eqs. (|T3p and (|T^). 

For /i^ ^ b^ (and /i^ ^ m^) the value of the mass function at zero momentum increases 
linearily with increasing renormalization point, 

M(0) oc /i , (33) 



a scaling behavior that is well known from studies of four-dimensional QED |20|. The vector 
part of the quark self energy. A, tends toward unity while the scalar part B admits an almost 
constant value of order // up to s ~ /i^/2. Obviously, the DSE-model defined by Eqs. (^ 
and dTDI ) is just an IR-regularised version of four- dimensional QED for b/^ <^ 1. 



Modifying the asymptotic behavior of the gluon exchange by taking into account the 
renormahzation group (RG) improved loop-corrections, cf. Eq. (|^), introduces a third mass 
scale: Aqcd- Hence, no scaling comparable to that of Eq. (Q) occurs anymore. Instead, the 
quark self energy functions A and B show the proper RG transformation properties. For 
several parameter pairs (6/ Aqcd, A) we have studied the large momentum behavior of the 
quark mass function M and have found a very good agreement with an asymptotic analysis 
based on the RG. This is illustrated in Fig. 2 where the calculated mass function for A = 30, 
6/Aqcd = 1, and /x/Aqcd = Auv/Aqcd = 2^^ is shown and an asymptotic tail of the form 



s 
with 



, d-l 

M{s) "" 




(34) 



c= 1.642 A^CD and rf = 0.729 (35) 

is fitted. The coupling constant A extracted from this UV tail coincides within 5% with 
the value inserted into the DSE {d = -^ ^ 0.760) demonstrating good numerical accuracy. 
(Note that d is the exponent of a logarithmic term.) The quark condensate in the model 
characterized by the above parameter set can be calculated from Eq. ( pS] ) together with 
Eq. (^5|); its value is 

(gg)M^=iGev2 - (-(0.70 ± 0.03)Aqcd)' ~ - ((190 ± 60) MeV)=^ . (36) 

We have used Aq^^ = 270 ± 80 MeV ^2[ . The same result is obtained calculating the 
trace of the quark propagator, Eq. (^Ol): (qq) ^2^iq^y2 ^ (— (0.67 ± 0.02)Aqcd) • 

In Fig. 3 the order parameter M(0) is shown as a function of the coupling constant A. 
For vanishing current mass, m = 0, we have studied the dependence on the renormalisation 
point fj, (here set equal to the cutoff Auv) choosing the IR mass scale b equal to Aqcd as 
suggested by Refs. P,0. The critical behavior of the order parameter, M(0) oc (A — Ac)'' 
of the solution in the chiral limit gets softened if the current quark mass is non-vanishing, 
i.e. we find the well-known behavior of an order parameter of DCSB. While increasing the 
regularization parameter Auv as well as the renormahzation point /i the critical coupling 
Ac decreases slightly. The extrapolation Ayy -^ C)0 yields a value of the critical coupling 
of Ac = 28.98 lb 0.02. As in case of A > Ac the vector part of the quark self energy. A, 
is quite close to unity, the Nambu-Goldstone phase of DCSB is dynamically favored. The 
critical exponent extracted from the data is (3 = 0.59 ± 0.01. This value is almost identical 
to /? = 0.57 ~ 0.58 Hawes et al. have derived in their study of the DSE system (^ ^2[) 
inserting several dressed quark-gluon vertices . 

The critical coupling strongly depends on the IR scale b: Small values of 6/Aqcd lead to 
small, large ones to large critical coupling constants as one might expect. For Np = (2) 
the IR scale that yields a critical coupling identical to the coupling suggested by the QCD 
asymptotics, Ac = fvr^ (fvr^), is &7Vf=o (2)/Aqcd ~ 0.26 ± 0.01 (0.28 ± 0.02). Therefore 
only for very small IR mass parameters b < 100 MeV DCSB takes place. 

In Fig. 4 the mass functions for different IR parameters b are shown for A = ||vr^. The 
parameter b is varied by two orders of magnitude (for &/Aqcd > 0.3 no DCSB occurs). It is 
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remarkable that neither the strength of DCSB as "measured" by the order parameter M(0) 
nor the momentum range that contributes dominantly to the mechanism of DCSB depend 
on the IR scale b introduced by the effective gluon exchange, i.e. the scales of DCSB and 
the gluon mass scale b are not or only loosely connected to each other if a realistic coupling 
constant A is used. 



IV. ANALYTIC CONTINUATION OF THE QUARK SELF ENERGY 

A. Mass Shell Condition and Constituent Quark Mass 

To study quark confinement we will analyze the quark propagator for timelike momenta. 
Therefore the Euclidean quark self energy has to be continued from spacelike to timelike 
momenta p -^ ip. Correspondingly, the variable s = p^ has to be continued to negative 
values, s < 0. 

Eqs. (|13]) and (0) can be understood as the defining equations for the vector and scalar 
parts of the quark self energy, A and 5, respectively. Therefore the analytic continuation 
can be performed by smoothly changing the external momenta s while keeping the loop 
momenta r spacelike. Thus one generates the quark self energy as a function of complex 
momenta by simply using complex s in Eqs. (|T3p and (|T^. Obviously, A*{s*) = A{s) and 
B*{s*) = B{s) and hence solving the DSE in the upper half of the complex plane is sufficient. 
Note that the functions A{r) and B{r) inserted into the right hand side of Eqs. ( [TB| ) and 



(p!^) are still those determined self-consistently on the spacelike momentum axis. 

First we generate the functions A and 5 on a small strip including the positive real axis. 
There we find that A and B change continuously. No singularities appear and therefore 
the analytical continuation can be performed patching open sets towards timelike momenta, 
s < 0. Problems arise caused by the poles of the angular integrals / dx Ia/b{x; r, s)E{x; r, s). 
Those are kinematical singularities that appear for sq > ^^.^^^.^ where s = SQexp{ii!}) with 
So = |s|. They are integrable as can be seen from the analytic expressions ( [A 3D and ([A4|) of 
the appendix, for instance, or from a Laurent series expansion of the kernels ( |TT| ) and (|T^. 
Nevertheless, the numerical results for the domain sq > ^^.^^a^,^ have been quite unstable 
against a change of numerical parameters. Hence we restrict our discussion on results for 
So < 2siiv^i?/2 o^ly- Both functions, A and B, are entire on this domain. 

The identification of a mass pole in the quark propagator fullfilling condition (^Tj) is 
only possible if 6 > \/2M{—b'^/2). This implies that in the case of DCSB our numerical 
method of identifying such a pole is only reliable for IR mass scale parameters b larger 
than the created dynamical mass M(0). In this restricted parameter space we find that 
the mass shell condition (^Tj) for the quark is fulfilled for both models, with and without 
UV-improvement. I.e. close to the momentum —sm quarks propagate like stable particles 
with constituent mass M{—sm)- For parameter values b < ^/2M(0) we are not able to make 
any statements about the existence of a quark mass pole. 

In Figs. 5a and 5b the mass function M{s) in the UV-improved model (|^) for 
{b/AqcB,^) = (0.27, Aqcd = 11^^) ^^^ (1)30), resp., are shown. The parameters used 
are such that the kinematical singularities are far away from the timelike momentum region 
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considered here. One recognizes that the mass condition ( plD is fullfilled and that the pole 
mass yJSM is very well approximated by the dynamical mass M(0). 

Our result of quark deconfinement (obtained for such values of the IR parameter b where 
our numerical results are stable) is in full agreement with the results of Hawes et al. They 
have investigated the scalar part of the quark Schwinger function, A5(T), analyzed its 
asymptotic behavior for large Euclidean times T, and extracted the mass Mas of a stable 
asymptotic fermion state through the relation |]^ 

^lim Ain[As(r)] = -Ma3. 

Mas is found to be (positive) finite and of the order of Aqcd- While we are using the 
rainbow approximation (|T0|) for the quark-gluon vertex Hawes et al. employ non-trivial 
vertices. These are constructed in a way that they obey the Ward-Takahashi identity, 
carry no kinematical singularities, go versus the free vertex for free fields, and ensure the 
propagators transforming properly under the Landau-Khalatnikov transformations [|23|,|l|. 
In addition, one of the vertices used in Ref. |0 leads to multiplicative renormalizability ||24|| . 
Our study confirms the result of unconfined quarks obtained for an IR vanishing quark- 
quark interaction . We may also deduce that in this model the form of the quark-gluon 
vertex is of minor importance (at least, as long as it is free of singularities). 

B. Analyticity and Wick rotation 

In order to study the analytic properties of the quark propagator we employ Cauchy's 
integral theorem in the following way. Suppose that the kernels 

/^(^5 s) = ^^2(^)^^2(^) /^ ^^^^(^5 r, s) E{x; r, s) , (37) 

B(r) /•! 

of the integrals J a and Jb (cf. Eqs. ([T5|) and ([T6|)) are analytic. Then the integral of /a(b) over 
a closed contour in the complex momentum plane has to vanish. Due to the renormalization 
condition ([T8|) the functions fA{B) fall off faster than 1/r^ for large r where x > 1. Therefore 
one obtains 

T"' drfAiB) (r) ^ T"" dr e' Va(b) (r e^') (39) 

where the difference between the left hand side and the right hand side vanishes for 
Auv —>■ cxo. This relation implies that the iteration of Eqs. (|T3|) and (|T^) along a rotated 
axis should lead to the same values A{0) and -B(O) as the iteration along the spacelike axis 
as long as the kernels fA{B) are analytic. This provides us with an analyticity test for the 
functions /a(b)- 

The gluon propagator (^ is symmetric under the exchange of spacelike and timelike 
momenta: 

Dri-k^) = -Drik^) . (40) 
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Even though the UV-improved version of this gluon exchange does not carry the symme- 
try (^0|) an UV-improvement with the replacement (H) does. As the spacehke properties 
(such as dynamical chiral symmetry breaking) are only slightly affected by the replace- 
ment (§) we believe our symmetry considerations being quite general and we concentrate on 
the UV-unimproved form (^ for the gluon propagator alone. 

In the following we study the approximate DSE system (^,^) and we drop the depen- 
dencies on regularization parameter Auv and renormalization point /i. Assuming J'a and Jb 
being analytic within the whole complex momentum plane the iterations of Eqs. (^) and 
(E^) along the timelike axis 



Gtt'^ Jo rA'^(—r) — B^(—r) J-i 



/l(-s) = l + — ( dr ^^^ ; L^^ ^ I dxlAix;r,s), (41) 



B{-s) = (m -^ 0) + A /^^^ dr ^^ "^ /' dx^x; r, .) (42) 

bn-^ Jo rA^{—r) — B^{—r) J-i 

should lead to the same values ^4(0) and -B(O) as the iteration along the spacelike axis. 
Obviously, 

A{-s) = A{s) , B{-s) = ±tB{s) (43) 

are solutions of Eqs. ( |H| ) and (^2|) if A and B are solutions of Eqs. (^Tj) and ([2^). In the 
Nambu-Goldstone realization of the vacuum, -B(O) 7^ 0, the solutions (|43|) do not obey the 
analyticity condition A{0) = A{0) and -B(O) = -B(O). Therefore the kernels /a and fs are 
not analytic in the whole complex momentum plane in case of chiral symmetry being broken 
dynamically. Note that a "usual" mass pole of the quark propagator is sufficient to explain 
the behavior (|^). 

The functions /^ and fs are products of the vector or scalar part of the quark propagator 
and the angular integral over Euclidean projections of the gluon propagator. We will argue in 
the following that the non- analyticity of both functions are caused by the quark propagator. 
For doing this we study the Wigner-Weyl phase of the system, B{s) = 0. In that case the 
self-consistency equation for A{s) = ai(so,'i9) + i 0:2(50, f?) reads: 

ai{so,i9)+ia2iso,i^) = 1 + —^ dr— — — -— — / rfx/^ (x; r, sqc*'') (44) 



with a\ and 02 being the real and the imaginary part oi Ais = sqc*''), respectively. It is 
straightforward to proof the following symmetries 

"1/2(5, TT + ^) = ai/2(s,^) , (45) 

ai(s,^-^) = ai(s,t9), (46) 

02(5,^-^) = -02(5,^). (47) 

Hence it is sufficient to analyze the integral equation ( PD for -Q G [0, |] only. 

In Fig. 6 the function a;i(so, "&) is shown for A = 16 and various angles d. We have used a 
small cut-off, A^y = 2^6^. However, the solutions do not depend on a further increasement 
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of AuY as the chiral symmetric solutions do not display the scaling behavior (|33|) as in the 
case of the Nambu-Goldstone phase. The numerical error in Fig. 6 is about 1% to 5%; it is 
largest close to the imaginary axis. 

It is obvious that for all the angles •& the real parts ai at zero momentum {s = 0) are 
identical; the imaginary parts 0:2 which are not shown vanish for small momenta within the 
numerical accuracy. For i) = ^ = 90° the iteration of Eq. ( ^If ) does not converge due to the 
(integrable) singularities of the gluon propagator. Nevertheless there is no signature of any 
further non-analyticity close to the imaginary axis. We therefore conclude that there is a 
strong indication for the vector part of the quark self energy being analytic. Then the quark 
propagator has only one pole, namely the one at s = 0. 

In the Nambu-Goldstone phase, B{s) 7^ 0, the integral kernels /a and fs from Eqs. (|37D 
and p8| ) are not holomorphic on the whole complex momentum plane either. There has to 
be a pole contribution from the contour integration in Eq. (|39|) . As the angular integrals 
over Ia and Ib do not carry a pole singularity the scalar and the vector parts of the quark 
propagator have at least one pole. If this singularity is not lying on the real axis but within 
the complex plane the Wick rotation cannot be performed naively. If it is placed on the 
timelike axis the naive Wick rotation is allowed as long as the causal behavior of the quark 
propagator is suitably determined. 

We have not found any singularities except the kinematical ones and the mass pole on 
the real axis. Even though we are not able to draw compulsory conclusions this numerical 
result strongly indicates that the quark propagator is analytic within the whole complex 
momentum plane except on the real timelike axis. This statement is also supported by 
the analysis within different DSE models [^ where it is found that poles of the quark 
propagators are lying quite close to the real axis. 

In addition, the analyticity condition ( P5| ) is quite easily understood from a physical 
point of view. Under the action of the Wick rotation the Euclidean mass term becomes (— z) 
times the corresponding Minkowski term, a kinetic fermionic term keeps "unmodified". In 
terms of the quark propagator functions this reads: 

B{sE)'^-tBisM), (48) 

A{se) ^' A{sm) (49) 

where se denotes the Euclidean and sm the corresponding Minkowski four momentum 
squared close to a mass pole of the quark propagator. In Minkowski space we therefore 
just recover one of the "solutions" ( P3D which is not surprising. 

Summarizing this section, there are strong indications that the Euclidean quark prop- 
agator of the DSE model specified by the gluon exchange (|^) can be continued back to 
Minkowski space. Within a small region around the mass pole (0) the Minkowskian quark 
propagator has a Kallen-Lehmann representation and it therefore resembles the propaga- 
tor of a massive stable fermion. For a restricted set of parameters we have been able to 
verify this numerically. Due to the numerical problems associated with the kinematical 
singularities we are not really able to prove this statement rigorously for all values of the 
parameter b. Nevertheless, it is the most plausible conclusion drawn from the properties of 
the functions f^ and fs as discussed in this section. 
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V. SUMMARY AND CONCLUSIONS 

We have investigated the Euchdean Dyson-Schwinger equation (DSE) for the quark self 
energy for an infrared vanishing gluon propagator in rainbow approximation. We have found 
that this model is not appropriate for describing hadron phenomenology due to several 
reasons. First, dynamical chiral symmetry breaking (DCSB) occurs either for unphysically 
large coupling constants only or for a small IR scale, b <^ Aqcd- The strength and the 
scale of DCSB are then independent of b: the spontaneous breakdown of chiral symmetry is 
mainly driven by the term added to obtain the correct UV behavior and not by the infrared 
part. 

In addition, continuation to timelike momenta reveals that the quark propagator pos- 
sesses a pole. As our numerical method is only reliable for IR mass scale parameters b larger 
than the created dynamical quark mass (see the discussion in subsection PVA| ) this pole can 



be unambigously identified for such values of b only. To be explicit, we do not claim that 
there is no pole, we only state that within our method it cannot be found even if it exists. 
However, the occurrence of a mass pole in the propagator can be related to the propaga- 
tion of a stable particle. Therefore it is clear that quarks are not confined in this model, at 
least for a large range of the IR parameter b. This is in contradiction to results obtained in 
Ref. |2^. There heavy quarkonia have been studied using a non-relativistic reduction of the 
gluon propagator (j^) together with a quark propagator and a quark-gluon vertex which re- 
sult from the DSE studies of Ref. [Q]. However, in this work the quark-quark interaction (^ 
has been replaced by a Coulomb potential. The resulting equation for the quarkonium 
bound states is just a Lippmann-Schwinger equation for hydrogen with the exception that 
the non-trivial quark propagators of Ref. [ffl have been used. These propagators are then 



probably the reason for the "confining mechanism" described in Ref. p5| . Stated otherwise, 
the unusual behavior of the quarks have been put in "by hand" . Such a mechanism is not 
supported by our calculations. 

Instead, our results are in accordance with those of Hawes et al. ^ who have studied 
the quark propagator for spacelike momenta only and have extracted the quark mass from 
the large time behavior of the Schwinger functions. Additionally, they have used non-trivial 
quark-gluon vertices. The good quantitative agreement of our results with those of Ref. ^ 
shows that the rainbow approximation is a good approximation for this model interaction. 
This can be understood from the fact that the influence of a non-trivial quark-gluon vertex 
is concentrated in the infrared which for the studied case is suppressed by the infrared 
vanishing gluon exchange. 

Even though a study of the analyticity of the quark propagator is hampered by numerical 
difficulties we have found convincing evidence for the quark propagator being analytic within 
the whole complex Euclidean momentum plane except the timelike (negative real) axis. This 
has two consequences: First, it is possible to perform the naive Wick rotation. Second, and 
more important, the model does not lead to quark confinement. 
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APPENDIX: 

The angular integrands in Eqs. (p!3D and ([T^) are: 



2 1 + 2x^)^/rs - 3a; r + s) 



Ux,r,s) ^ \I-Vl^- \\ '\ 'I ' , (Al) 

[r + s — Ix^rs) + 0* 

3(r + s — 2xJTs 



Mx;r,s)^rVT^^ \ ' ^ ^. ' . (A2) 

(r + s — Ix^rs) + r 



For S(a;; r, s) = 1 in Eqs. (|13[) and (|1^ {i.e. no UV improvement for the quark-quark 
interaction) the x-integral can be evaluated exactly: 



/ dxlA{x]r,s) 



Td^ _ 7nf_ r 2 ( jr-sf + b^ _ (r + g)^ - b^ {r + s)\{r - sf + b^) \ 
4^ ~sy' { V(r + s)2 + 64 2^s ^ rs{{r + s)"^ + b^) / ^^ 

(r + .)((r-.)^ + 5&^) 1 

H 7 T^ — 71 ^- } 5 (A3) 

(r + s)2 + 64 J ' ^ ^ 

l\dxlnix;r,s) = ^^^^^^ (A4) 



((r + s)2 + b^)y I Ars 
with 

y = {{r- sf + by + IQrsb^ + ((r - sf + b^)J{{r-sy + by + 16rsb^, (A5) 



w± = \/ ((r + s)2 + b^)^{{r - s)^ + b^ + IGrsb"^ ± (((r + s^ + b^ - 4rs((r + s)^ - b^)). 

(A6) 



16 



REFERENCES 

EMaii alkofer@ptsunl.tphys.physik.uni-tuebingen.de, 
bender@theory.phy.anl.gov 
[1] CD. Roberts and A.G. Williams, Dyson-Schwinger equations and their Application to 
Hadronic Physics^ Progress in Particle and Nuclear Physics, Vol. 33, pp. 477-575, ed. 
A. FaBler, Pergamon Press, Oxford, 1994. 
[2] Baryons are described by relativistic Fadde'ev equations. 
[3] D. Zwanziger, Nucl. Phys. B 364, 127 (1991). 

[4] U. Habel et al, Z. Phys. A - Atomic Nuclei 336, 423 and 435 (1990). 
M. Stingl, Phys. Rev. D 34, 3863 (1986); Erratum zbid. D 36, 651 (1987). 
M. Stingl, A Systematic Extended Iterative Solution for Quantum Chromdynamics, 
Universitat Miinster preprint MS-TPI-94-13. 
[5] L. V. Smekal, Ph.D. thesis, Tiibingen University, 1995. 
[6] C. Bernard, C. Parrinello, and A. Soni, Phys. Rev. D 49, 1585 (1994). 
[7] F.T. Hawes, CD. Roberts, and A. G. Williams, Phys. Rev. D 49, 4683 (1994). 
[8] We use negative semidefinite Euclidean metric, gfj^^y = —S^u, a ■ b = — Z]j=i cijbj < 
Wa = b E IR^. For abbreviation, we nevertheless define the square of a momentum 
as q^ := —q ■ q > 0. All Euclidean Dirac matrices, 7^, are antihermitean and obey the 
anticommutation relations {7^,71.} = —26^^,. 
[9] W. Marciano and H. Pagels, Phys. Rep. 36, 137 (1978). 
[10] L. P. Pascual and R. Tarrach, QCD: Renormalization for the Practitioner, Lecture 

Notes in Physics 194, Springer- Verlag, Berlin - Heidelberg - New York, 1984. 
[11] R. Delbourgo and M. D. Scadron, J. Phys. G 5, 1631 (1979). 
[12] C D. Roberts, private communication. 

[13] R. T. Cahill and C D. Roberts, Phys. Rev. D 32, 2419 (1985). 
[14] H. D. Politzer, Nucl. Phys. B 117, 397 (1976). 
[15] F. J. Yndurain, The Theory of Quark and Gluon Interactions, Springer- Verlag, Berlin 

- Heidelberg - New York, 1993. 
[16] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields, McGraw-Hill, New York, 

1965. 
[17] J. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D 10, 2428 (1974). 
[18] C D. Roberts and R. T. Cahill, Phys. Rev. D 33 1755 (1986). 
[19] D. Atkinson and P. W. Johnson, Phys. Rev. D 35 1943 (1987). 
[20] R. Fukuda and T. Kugo, Nucl. Phys. B117, 250 (1976). 

[21] S. J. Stainsby and R. T. Cahill, Phys. Lett. A 146, 467 (1990); Int. J. Mod. Phys. A 7 
7541 (1992). 

P. Maris and H. Hohies, Int. J. Mod. Phys. A 7, 5369 (1992). 
[22] Particle Data Group, Phys. Rev. D 50, 1173 (1994). 
[23] L. D. Landau and I. M. Khalatnikov, Sov. Phys. JETP 2, 69 (1956). 
[24] D. C Curtis and M. R. Pennington, Phys. Rev. D 42, 4165 (1990). 
[25] M. Becker et al, Phys. Lett. B 267, 261 (1991). 



17 



FIGURES 

FIG. 1. The scalar part of the quark self energy, B{s)/B{0), is shown as a function oi s = p^. 
The UV-unimproved gluon propagator (^), the coupling strength A = 16, and the numerical cut-off 
Ayy = 2^^b^ have been used. The solid line represents the function obtained by iterating Eqs. (13) 
and (jlj). The dashed one is the asymptotic tail (26) fitted at s = 2^^ 6^ with the parameters (32) 
inserted. 

FIG. 2. This figure shows the quark mass function M(s) of the DSE-model with UV-improved 
gluon propagator (0). The parameters used are: A = 30 and 6/Aqcd = 1- The solid line represents 
the function obtained by iterating Eqs. (^) and (p!^. The dashed one is the asymptotic tail (34) 
fitted at large momenta s with the parameters ( |35|) inserted. 



FIG. 3. In this figure the order parameter M(0) is shown as a function of the coupling strength A 
for the model with the UV-improved gluon propagator (Q) and 6/Aqcd = 1- The thick lines 
represent solutions of the DSE (13) and (^) for vanishing quark current mass, m = 0, and different 
renormalization points /i^ = Auv- The thin lines are solutions for different current masses and 
fixed /iVA^jcD = A?jv/A^cD = 2^'- 

FIG. 4. In this figure the quark mass functions for different IR parameters b are shown. The 
UV-improved gluon propagator is used and the coupling strength is A = ^'^'^- The momentum 
range that contributes dominantly to the breaking of chiral symmetry is defined as the value of s 
where M{s) = ^M(O). 

FIG. 5. The quark mass function at spacelike and timelike momenta 
for (6/Aqcd,A) = (0.27, Aqcd = 48/(297r^)) and (1,30) are shown in Figs, a and b, resp. The 
thin line on the timelike side of the graphs is the function \/—s. The point of intersection of this 
function with the mass functions gives the value of s where the pole condition (pTl) is fulfilled. 



FIG. 6. In this figure the real part ai of the scalar quark self energy A is shown for the 
Wigner-Weyl phase of the DSE model defined by the gluon exchange (P). The different functions 



correspond to solutions of Eq. (44) iterated along different axes in the Euclidean momentum plane 
characterized by the angle ^. 
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